Mon. Not. R. Astron. Soc. 000, [l[|| (1994) Printed 5 February 2008 (MN WT$L style file vl.4) 



Magnetically modulated accretion in T Tauri stars 

C.J. Clarke*, P.J. Armitage, K.W. Smith and J.E. Pringle 

Institute of Astronomy, Madingley Road, Cambridge CB3 OH A 



Accepted 1994 October 26. Received 1994 August 18 



in 

On 
(X, 



O 
(N 
O 

in 

Oh 
i 

o 



13 



ABSTRACT 

We examine how accretion on to T Tauri stars may be modulated by a time-dependent 
'magnetic gate' where the inner edge of the accretion disc is disrupted by a varying 
stellar field. We show that magnetic field variations on time-scales <§; 10 5 yr can 
modulate the accretion flow, thus providing a possible mechanism both for the marked 
photometric variability of T Tauri stars and for the possible conversion of T Tauri stars 
between classical and weak line status. We thus suggest that archival data relating to 
the spectrophotometric variability of T Tauri stars may provide an indirect record of 
magnetic activity cycles in low-mass pre-main-sequence stars. 
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1 INTRODUCTION 

The relationship between classical and weak line T Tauri 
stars (henceforth CTTs and WTTs) remains an open ques- 
tion in the study of pre-main-sequence stellar evolution. 
Whilst a number of spectral diagnostics suggest that CTTs 
are associated with material close to the stellar photosphere, 
it is not all clear whether the weak line systems (which lack 
such signatures) are necessarily devoid of circumstellar ma- 
terial at all radii and are thus, permanently and irreversibly, 
truly naked systems. The strongest argument that this is not 
the case is provided by the observation that a number of T 
Tauri systems have been observed to switch betwe en strong 



and weak line status on a time-scale of decades (Herbig & 



Bell 198c). Since this time-scale is considerably less than 



even the fastest (dynamical) time-scale associated with any 
reservoir of material at large radius it is clear that, in these 
systems, the transition between CTT and WTT status is 
controlled by processes close to the star. By implication, 
therefore, at least some stars classified as WTTs must be 
indistinguishable on the large scale from some CTTs. Fur- 
ther support for this view is provided by the overlap in stel- 
lar parameters between the WTTs and CTTs as a class, 
both in terms of their location on the Hertzprung-Russell 
(HR) diagram and in terms of their angular momentum per 
unit mass. Furthermore the detection of spectral excesses in 
some WTTs at submillimetre or near-infrared wavelengths 
has been claimed as evidence that some systems are 'weak' in 
terms of their circumstellar material close to the star but are 



nevertheless associated with matter at large radii (Strom et 
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al. 1989; Beckwith ct al. 199C|), a lthough this claim remains 
controversial ( |Waltcr et al. 1992[ ). 

In this paper we propose a simple mechanism that can 
allow systems to switch between weak line and strong line 
status, and may also give rise to the pronounced photomet- 
ric variability of T Tauri stars on a time-scale of decades. We 
argue that variations in both the lines and continuum are re- 
lated to accretion processes close to the star, and that these 
are modulated by a time-dependent magnetic 'gate' at the 
inner edge of an accretion disc. In this model the disc is dis- 
rupted by the magnetic field of the central star at a distance 
of a few stellar radii; material can, however, only flow down 
the field on to the star if the star is rotating more slowly 
than the local Keplerian frequency of the disc at the mag- 
netosphere. If, therefore, the field is variable, the changing 
location of the magnetosphere can give rise to alternating 
phases of storage of material in the disc and of accretion on 
to the star. We thus suggest that the variability pattern of 
T Tauri stars may provide indirect evidence of the cycles of 
magnetic activity in pre-main-sequence stars, analogous to 
the similar record infer red from chrom ospheric variations of 
main-sequence dwarfs (Baliunas 1988). 

In Section 2 we set out the model in more detail, es- 
timating the field strengths and variability time-scales for 
which the magnetic gate would work in the way outlined 
above. In Section 3 we discuss the model in relation to a 
variety of observational parameters of T Tauri stars. 



2 THE MODEL 

In this section we set out the model of magnetically coupled 
T Tauri-disc systems and derive the conditions under which 
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such a model can modulate the accretion flow on to the diffusion in the disc (Campbell 1992; flomanova, Lovelace 



central star 



Bisnovatyi-Kogan 1995) or else the action of a poorly 



We assume that the central star possesses an ordered 
magnetic field on large scales whose axis is aligned with the 
rotatio n axis of the star. Furthermore, following Ghosh & 
Lamb (1979) and subsequent authors we assume that the 



field is able to penetrate the disc as a result of a variety of 
micro-instabilities that ensure efficient mixing between the 
field and disc plasma. Apart from a narrow region around 
corotation between disc and star, the mismatched angular 
velocity between field lines (corotating with the star) and 
Keplerian disc material ensures that the field exerts a torque 
on the disc, transferring angular momentum from/to the disc 
respectively inward/outward of corotation. 

In such a situation, a disc of material flowing inwards 
under the action of viscosity is first significantly perturbed 
by the field at radius r m where the magneti c torque becomes 



comparable wit h the local viscous torque (Bath et al. 1974 



Campbell 1992). (Note that, owing to the ordering of the 
magnitude of the viscous torque per unit volume and the 
thermal and dynamical energy densities in accretion discs, 
such a criterion is satisfied when the magnetic pressure is 
negligible compared with either the thermal or ram pressure 
in the disc; one can therefore use unperturbed disc structure 
equations in the estimation of r m .) For steady-state disc 
flow around a star of mass M„, at accretion rate M, the 
differential torque across an annulus of width dr is given by 



dG„ 
dr 



M ( GM* 
~~2 \ r 



1/2 



(1) 



where G„ denotes the torque due to viscous processes in the 
disc. The corresponding differential magnetic torque is 

dGr, 



dr 



rB^B 



(2) 



where B z and B$ are the vertically averaged poloidal and 
toroidal components of the field respectively. Equating these 
quantities, the magnetospheric radius then satisfies 



BzBjy 



M ( GM, x ' 1 



(3) 



We follow previous authors (e.g. Cameron & Campbell 1993; 
Konigl 1991) in assuming that, whatever the chaotic field 
structure near the stellar surface, the residual poloidal field 
at several stellar radii is approximately dipolar: 



B z 



B, 



(4) 



(where B* is the equatorial dipole field at the stellar sur- 
face, r„) and that, moreover, the equilibrium poloidal field 



threadi ng the disc is close to it s unperturbed value (Ghosh 
& Lam fj 1979| ; |Campbell 1992] ). Evaluation of the toroidal 
field, however, requires several assumptions about the pro- 
cesses determining the equilibrium value of B^ in a shearing 
velocity field. Whereas the mechanism for the generation 
of toroidal field is uncontroversial (i.e. differential rotation 
between the footpoints of field lines anchored respectively 
on the star and disc), a number of mechanisms have been 
suggested for limiting the growth of B^, including r econ- 
nection of twis ted field lines in the magnetosphere (Livio 
& Prinjde 1992), magnetic buoyancy or turbulent magnetic 



specified anomalous resistivity, whose role is to dissipate B^ 
on a time-scale comparable to the A lfven crossing time of 
the disc height (Ghosh & Lamb 1979). Here we assume that 
magnetic reconnection in the magnetosphere is the domi- 
nant limiting process, which implies B$ ~ B z since, away 
from corotation, B$ is both generated and destroyed on a 
dynamical time-scale. Thus equation (H) becomes 



B, L = -— 



M ( CM. N 1 1 



(•>) 



It is immediately apparent that this expression is within a 
numerical factor of order unity of that obtained if one calcu- 
lates the field strength required to dominate magnetically a 
spherical inflow of the same M (Davidson & Ostriker 1973), 
thus justifying, to order of magnitude, earlier works that 
have applied the spherical Alfven radius in the context of 
disc accretion (e.g. Konigl 1991). A particularly convenient 
aspect of this result is that r m depends only on the field 
strength and M: 



\{GM*y/ 2 M 



2/7 



(6) 



and is independent of the internal structure of the disc giving 
rise to this steady-state accretion rate M. 

We now consider the response of such a magnetically 
coupled star-disc system to cyclical variations in the field 
strength of the central star. According to equation ^ such 
changes produce cyclical variations in the magnetospheric 
radius r m . At strong field phases of the cycle, r m may lie 
outside the radius of corotation between disc and star, rn, 
so that the action of the field is then to impart angular 
momentum to the disc. At such phases, therefore, material 
is held up in the disc by the strong magnetic torques at its 
inner edge and is unable to flow down the field on to the 
star. As the field subsequently reduces, material can again 
flow inwards under the action of viscosity; once inside vq the 
magnetic torques change sign so that disc material can be 
decelerated and again flow on to the star. Thus modulation 
of the field can give rise to an intermittent accretion pattern, 
with strong field phases in which material is stored in the 
disc alternating with weak field phases in which the stored 
material flows on to the star. 

Several criteria have to be satisfied, however, in order 
for this process to operate in the way outlined above. Since 
the modulation of accretion depends on the alternate dis- 
placement of r m inward and outward of corotation it is nec- 
essary, for order-unity field variations to achieve this effect, 
that ra is close to the mean value of r m . Secondly, it is also 
necessary that the field variations are sufficiently rapid for 
variations of r m to occur at approximately constant rn. Be- 
low we consider each of these issues in the context of T Tauri 
stars. 

For a star rotating with period P* , the corotation radius 

is 

3 2\ 1/3 



GM,P j 

4tT 2 



(7) 



Thus (comparing equations ^ and it is evident that, for 
parameters typical of classical T Tauri stars (M ~ 10 -7 M* 
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yr -1 , P» ~ 8 d), rn and r m are of similar magnitude for 
surface field strengths of a few hundred gauss. Such field 
strengt hs are within the upper lim i ts derived from Zeeman 
effect ( Johnstone & Penston 1987; Ba sri, Marcy fc Valenti 
1992T 
1981 



and polarization measurements (Brown & Landstreet 



in T Tauri stars, thus giving support to the idea that 



T Tauri rotation is controlled by magnetically coupled disc 
braking. 

We now consider the maximum time-scale for field vari- 
ations that could give rise to intermittent material flow on 
to the star: if the field were to change sufficiently slowly 
then the star's rotation could re-adjust so as to maintain 
a fixed relation between rn and r m and hence maintain a 
steady flow from disc to star. For the case r m > rn (for 
which there is no accretion on to the star) the time-scale for 
stellar spindown, in, may be estimated as follows: assuming 
that field-disc interaction is concentrated over a region of ra- 
dial extent r m then (from equations S and the spindown 
torque exerted by the field is ~ Afr m f2 m , where £7 m is the 
Keplerian frequency at the magnetosphere. Thus 



ta 



fk 2 M t 
M 



(8) 



where k ~ 0.4 (van den Heuvel 1968) is the dimensionless 
radius of gyration of the star and /(< 1) is the ratio of the 
specific angular momentum of the stellar surface to that of 
material in Keplerian rotation at the magnetosphere. (Note 
that, for variations on time-scales <C 10 6 yr, the partial offset 
of this spindown by spinup due to the star's contraction on 
its Kelvin-Helmholtz time-scale may be neglected.) For a T 
Tauri star of rotation period ~ 8 d, r*/rn ~ 0.2 and thus, 
for r m modestly beyond rn, 



/ 



2 / N 1/2 

\ I rn 

rn 



0.04. 



(9) 



Hence, for a T Tauri disc with accretion rate ~ 10 -7 M* 
yr - , in ~ 10 5 yr. A similar argument may be made for 
the case where r m < rn, modified by the fact that in this 
case material flows from the disc on to the star. This in- 
troduces an additional angular momentum input rate of 
Mr^D» (which is, however, smaller than the spinup torque 
from the field by a factor f2*/n m < 1) plus an additional 
spinup torque due to the star's (adiabatic) contraction in re- 
sponse to accretion; this latter effect occurs on a time-scale 
~ M,/M ~ 10 7 yr. The dominant contribution to changing 
the stellar period is thus the magnetic torque (whether r m 
is greater or less than rn) and the appropriate time-scale 
is ~ 10 5 yr in either case. We therefore deduce that field 
variations on time-scales substantially less than 10 5 yr oc- 
cur at approximately constant rotation, and thus result in 
the modulation of the relative values of r m and rn. 



3 DISCUSSION 

In the previous section we showed that, for field strengths 
of a few hundred gauss, order-unity variations of the field 
on time-scales <C 10 5 yr would strongly modulate the ac- 
cretion flow on to the central star. Further characterization 
of the time-dependent flow and associated photometric vari- 
ability requires detailed hydrodynamic modeling (Armitage 
et al., in prep). The main implications of this conclusion are, 



first, that spectrophotometric variations of T Tauri stars 
may thus provide a record of the magnetic activity cycle of 
T Tauri stars and, secondly, that at least some weak line 
systems are interconvertible with strong line systems. We 
examine each of these issues below. 

The temporal behaviour of magnetic fields in T Tauri 
stars is not known, apart from the rapid variations deduced 
(by analogy with solar chro mospheric /coronal flares) from 
short-time-scale ultraviolet (Worden et al. 1981) and X-ray 
(Montmerle et al. 1983, 1993) flares in these objects. Further 
analogy with magnetic activity in main-sequence stars leads 
one to anticipate that T Tauri magnetic fields would also 
be variable on a range of longer time-scales, as witnessed by 
the chromospheric variations of late- type dwarfs (time-scales 
~ 10 yr: Baliunas 1988) and the history of intermittent ac- 
tivity on longer time-scales contained in the sunspot record. 
The qualitative similarity between the convective dynamos 
believed to be responsible for magnetic fields in both main- 
sequence and pre-main-sequence low-mass stars leads one to 
anticipate a similarly variable pattern in T Tauri stars. 

Pronounced variability, on a range of time-scales, is one 
of the defining characteristics of T Tauri stars as a class. 
Apart from the rapid flares alluded to above, which may be 
interpreted as localized surface phenomena on the star, it is 
hard to avoid the conclusion that this variability is related to 
variable accretion on to the star since, in the absence of ac- 
cretion, it is difficult to envisage what processes in the stellar 
interior could account for such large- amplitude variations in 
the bolometric output of the star on these time-scales. Vari- 
able accretion can result either from magnetic variations (as 
described here) or else from variable flow in the accretion 
disc, although it is not at all obvious why the accretion flow 
should be intrinsically intermittent on these time-scales. We 
note that the hypothesis that T Tauri variability reflects 
the stars' magnetic activity opens up the possibility of us- 
ing the photometric variability record of T Tauri stars as a 
diagnostic of magnetic activity cycles in T Tauri stars, in 
much the same way as chromospheric variations are used 
in main-sequence stars. However, whereas the latter record 
only extends back a few decades, and for a limited sample 
of stars that have been deliberately monitored for this pur- 
pose, published photometry for T Tauri stars, together with 
archival plate material, means that a database extending 
back over a century is available for a large sample of stars. 

We next turn to the question of the interconvertibility 
of strong and weak line T Tauri stars. The strongest evi- 
dence for the possibility of switching between the two states 
is pr ovided by the half a dozen systems in the Herbig & 
Bell 



1988 



catalogue, for which it is noted that their cur- 
rent designation as weak or strong line does not accord with 
the results of objective prism surveys in previous decades. 
Since these remarks do not relate to systematic monitor- 
ing programmes they cannot be used to derive duty cycles, 
but just illustrate that such transitions can occur on time- 
scales of decades. Further plausibility arguments about the 
possible interconversion of WTTs and CTTs depend on the 
overlap in properties (apart from those indicative of mate- 
rial flow at small radii) between the two classes of object. 
It is well known, for example, that CTTs and WTTs co- 
exist in the same region of the HR diagram (apart from 
a WTT-only region near the main-sequence), thus arguing 
against a simple picture in which a universal clock controls 
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an irreversible transition of stars from WTT to CTT status 
The rotational properties of CTTs and WTT s also overlap 



whe reas the WTTs rotat e faster as a class ( Bouvier et a^ 
1993; E|d wards et al. 1993 ), about half of them have angular 



momenta that overlap with the CTTs. Since the slow rota- 
tion of T Tauri stars is commonly attributed to a history in 
which they have been braked by magnetic coupling to asso- 



ciated discs (Konigl 1991; Cameron & Campbell 1993), the 



low angular momentum WTTs have probably undergone a 
similar history of disc braking. Such data suggest a picture 
in which objects classified as WTTs are a mixture of gen- 
uinely naked systems (devoid of material at all radii) and 
those that, retaining a reservoir of circumstellar material, 
are capable of being resuscitated as strong line systems. 



The above picture suggests some correlations between 
the diagnostics of material at small radii and the rotational 
properties of T Tauri stars. For example, if the low angular 
momentum T Tauri stars have evolved to a state of rota- 
tional quasi-equilibrium (i.e. if the time-averaged magneto- 
spheric radius is close to corotation, so that the the spinup 
and spindown torques exerted by disc material nearly can- 
cel) then the magnetic field should disrupt the disc out to 
larger radii in more slowly rotating systems, thus increasing 
the volume of the region of quasi-spherical infall. Since Ha 
emission is believed to arise from material flowing nearly 
radially near the star (e.g. Calvet & Hartmann 1992) one 
would therefore anticipate a larger Ha flux in slowly rotating 
systems of given accretion rate, since the solid angle inter- 
cepted by inflowing material is larger than in more rapidly 
rotating systems. On the other hand, the range of rotation 
periods among CTTs may indicate that these systems have 
not all evolved to a state in which the magnetosphere lies 
at corotation; in this case the time-averaged magnetosphere 
would lie outside corotation in more rapidly rotating sys- 
tems. Since (see Section 2) accretion on to the star does 
not occur when the magnetosphere is outside corotation it 
follows that, in a magnetic cycle, the fraction of the cycle 
for which accretion can occur is lower for rapidly rotating 
systems. Thus in this case, also, one would expect weaker 
diagnostics of accretion close to the star for more rapidly 
rotating stars. It is notable that these predictions are borne 
out by the observed decline in the ratio of Ha (indicative of 
material close to the star) to excess K magnitude (indicative 
of the reservoir of material at large r radii) with increasing 
stellar rotation (Bouvier et al. 1993). 
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Finally we note that such a model is not unique in pre- 
dicting a variety of spectral energy distributions according 
to the inner radius to which the disc extends. The possibil- 
ity of systems that are 'weak' on small scales and 'strong' 
on large scales has been recognized by various authors who 
have sought to identify spectral 'gaps' in some T Tauri stars 
as clearing of the inner disc by planetary formation (e.g. 
Marsh & Mahoney 1992). Such a picture does not, however, 
unlike the variable magnetic disc model described here, al- 
low the reversible exhaustion/replenishment of material on 
small scales as suggested by the variability data described 
above. 
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